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Abstract 

The use of objectively analyzed fields of meteorological data for the 
initialization of numerical prediction models and for complex diagnostic 
studies places the requirements upon the objective method that derivatives of 
the gridded fields be accurate and free from interpolation error. We propose a 
modification of an objective analysis developed by Barnes that provides 
improvements in analyses of both the field and its derivatives. Theoretical 
comparisons, comparisons between analyses of analytical monochromatic waves, 
and comparisons between analyses of actual weather data are used to show the 
potential of the new method. The new method restores more of the amplitudes of 
desired wavelengths while simultaneously filtering more of the amplitudes of 
undesired wavelengths. These results also hold for the first and second 
derivatives calculated from the gridded fields. Greatest improvements were for 
the Laplacians of the height field; the new method reduced the variance of 
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undesirable very short wavelengths by 72 percent. Other improvements were 
found in the divergence of the gridded wind field and near the boundaries of 
the field of data. 

1. Introduction 

In the analysis of weather data, it is often necessary to transform the 
data from nonuniformly distributed observation sites to a rectangular grid 
suitable for numerical calculations. The transformation is accomplished by 
some form of interpolation and is referred to as "objective analysis". An 
important class of objective analysis techniques known as "successive 
corrections methods (SC)" originated into meteorology in 1955 (Bergthorsson 
and Doos, 1955) and was introduced into operational meteorology by Cressman 
(1959). Though recently replaced by multivariate statistical techniques as the 
operational interpolation method in some synoptic scale numerical forecasting 
models, the successive corrections methods remain among the most commonly used 
in practice (Bengtsson et al. , 1981). 

Because the SC methods have been widely used for gridding meteorological 
data, it is appropriate that they should be benchmarks for comparisons with 
newer objective analysis techniques (Seaman, 1983; Schlatter, 1976; Davidson, 
1982). It is therefore important that the intercomparisons be done with SC 
methods for which the analysis parameters have been set to yield as closely as 
possible the optimum analysis. Further, SC methods have been shown to be 
interchangeable with statistical interpolation methods (Bratseth, 1986; 
Lorenc, 1986), and since the SC methods are more economical, they will find 
greater application for the initialization of numerical prediction models. And 
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because the numerical models (and also some diagnostic, models) involve higher 
order calculations and products therefrom, it is crucial to obtain accurate 
derivatives from objectively gridded fields of meteorological data. 

This study presents a modification of the Barnes (196A, 1973) SC method 
to increase the accuracy of the objective analysis of smaller scale features 
contained in the data and to increase the accuracy of the derivatives 
calculated from the gridded fields. The Barnes method was chosen for this 
study because it has been widely used for gridding meteorological data and has 
been made part of data assimilation and analysis packages (Koch, etal. , 1983; 
Smith and Leslie, 1984). Some of the advantages of the scheme for this study 
are: (1) Its Gaussian weighting scheme has been incorporated into statistical 
methods (Seaman, 1983), (2) the weight function is not arbitrarily truncated 
(a practice that introduces high frequency noise into the gridded fields), (3) 
the response characteristics of the method have been well-documented, and (4) 
the shape of the weight function is not tied to an influence radius; the 
influence of data can be extended without changing the response 
characteristics of the weight function. 

Achtemeier (1987) results suggested that the desired improvements were 
possible if the traditionally used 2-pass version developed by Barnes (1973) 
were modified to include a third pass. Since the new method incorporates much 
of the original methodology developed by Barnes, we adopt the terminology, 
Barnes Objective ANalysis Three Pass (B0BAN3) for the 3-pass version. In the 
next section, we derive the response equations for B0BAN3 and demonstrate the 
theoretical increase in accuracy over the 2-pass method. In addition, 
comparative results from analyses of monochromatic sine waves are presented. 
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In Section 3, analyses of a jet streak embedded within a intense synoptic 
scale cyclone on 12 GMT 11 April 1979 are presented to show how B0BAN3 
improves over the 2-pass method. The significance of these results is 
discussed in Section 4. 

2. Theoretical Analysis 

Only portions of the theory developed by Barnes are reproduced here to 
provide the reader with the necessary background to understand this study. A 
more complete explanation is available in Barnes (1973) and is also contained 
in Koch etal. (1983) and Achtemeier (1987). 

The horizontal distribution of an atmospheric variable can be described 
mathematically by a linear combination of an infinite number of sinusoidal 
waves, a Fourier Integral representation. If we wish to suppress some of the 
information in the actual field and represent only components larger than an 
arbitrarily chosen wavelength, we can apply a judiciously-selected filtering 
function to the distribution. The effect upon the resulting depiction of the 
field can be demonstrated for a single sinusoidal wave of the form f(x,y) = A 
sin(ax), where a = 2 tt/L and L is the wavelength. The filtered (or objectively 
analyzed) field, g(x,y), is represented by 

n ® 

E(x+rcos9,y+rsin0)w(r ,k)rdrd0. (1 ) 

where the filter is a function of the distance, r, between the arbitrary point 
(x,y) and all other points in the field. Barnes (1964) chose a Gaussian 
function for w(r,k) 



(5) 


w(r,k) = (4Trk) - ^exp(-r^/4k) . (2) 

If the data concerning the atmospheric variable are distributed continuously 
over an infinite plane surrounding (x,y), as they are in this analytical 
example, then (1) becomes 

g(x,y) = D(a,k) f(x,y), (3) 

where the amplitude response, D(a,k), is given by 

D(a,k) = exp(-a^k). (4) 

Here k is a selectable parameter that determines the shape of the weighting 
curve for the filter. 

Equation (4) gives the response of the Gaussian filter after one 
application. Barnes (1964) showed that more of the details inherent in the 
data could be recovered if the filtered residual differences between f(x,y) 
and g(x,y) were added to the original filtered field. However, several 
additions of smoothed corrections to a smoothed initial analysis, were 
required to retrieve the important details inherent in the data. In 1973, he 
modified his method so that most of the important details are recovered in 
only two passes through the data. If k = kg and D(a,kQ) = Dg for the initial 
pass, and k^ = tkg for the correction pass, then the final response D' is 
given by 
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D' = D 0 (l + D 0 T_1 - D 0 T ). (5) 

Here x is an arbitrary weighting parameter that decreases the number of 
iterations required to achieve desired response. 

Achtemeier (1987) found from theoretical studies that equivalent or 
better accuracy could be obtained over the 2-pass 1973 (traditional) method by 
a 2-pass method with the same filter parameters on the initial and correction 
passes (fixed method). These studies required the final responses for the 
traditional and fixed methods to be equal at an arbitrary "reference 
wavelength" (usually the 2S wavelength where S is the average separation 
between observations). Accuracy was defined by how well the method retains 
desired wavelengths and removes undesired wavelengths. Applications revealed 
that the more sparsely and irregularly distributed are the data, the less the 
results are in accord with the predictions of theory. No single set of filter 
parameters for the 2-pass method produce the most accurate objective analysis 
for all wavelengths. The traditional method developed the most accurate 
objective analysis of the longer wavelengths. The fixed method added 
information present in the shorter wavelengths without increasing the 
contribution by undesired very short wavelengths. 

This information about how different filter parameters modify the 
performance of the 2-pass objective analysis was used to develop the 3-pass 
method that includes the desirable properties of both the traditional and 
fixed methods. This modification is but a minor extension of the existing 
technique. A very smooth first pass captures the long wavelengths and an 
additional 2-passes with the fixed method retains the desired shorter 
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wavelengths. Consider a 3-pass method for which the first pass response is Dq 
and the second and third pass responses are Dp The final response after 
3-passes is 

D' = D 0 + (1 - Dq^U + D 1 t_1 - D^) (6) 

where the third pass response is related to the second pass response by D 2 = 
rD^ . If Dq = 0, the first pass only restores the mean of the data field and 
has zero response for all wavelengths. Then the final response reduces to the 
response for the traditional 2-pass method given by (5). The final response 
for B0BAN3 is obtained upon setting x=l, 

D' = 1 - (1 - Dq)( 1 - D x ) 2 (7) 

Choosing a small Dq (say Dq = 2.5 x 10"^) allows the restoration of some of 
the longer wavelengths but has near zero response at the reference wavelength. 
Then the filter parameters for the fixed method give the desired final 
response at the reference wavelength after two correction passes. Fig. 1 shows 
the theoretical final response curves for B0BAN3 and the traditional method. 
The responses are identical (0.25) at the reference wavelength (2S). B0BAN3 
offers an additional 7-9 percent restoration of the 4S - 5S wavelengths. 

Though Dq has near zero response at the reference wavelength, and 
therefore contributes negligibly to the final response there, its contribution 
to the final response at the longer wavelengths is not negligible. The solid 
line in Fig. 2 shows the differences between the final responses for the fixed 
method and B0BAN3 expressed as percent of the amplitude of the original wave. 



( 8 ) 


The standard of comparison is the fixed method. The negative percentages mean 
that inclusion of the third pass causes restoration of more of the original 
wave. There is no difference between the two methods at the reference 
wavelength (2S wave). B0BAN3 restores two percent more of the 4S wave. It 
restores lesser percentages at the longer wavelengths because the 2-pass fixed 
method has almost perfectly restored these wavelengths already and there 
remains almost no additional information to add to the fixed analysis. 

Following a procedure for the analysis of monochromatic sine and cosine 
waves with analytical data (Achtemeier, 1987), we performed objective analyses 
by B0BAN3, the fixed method, and the traditional method for waves that range 
from 2S to 12S in wavelength. RMS differences between the true field and each 
analysis were calculated for the grid interior. The statistic of relative 
performance is the difference between the RMS differences for B0BAN3 and the 
traditional method and for the fixed method and the traditional method 
normalized by the amplitude of the analytical wave. It is a measure of how 
well each method restores the whole wave, not just the amplitude. This 
statistic approximates the percent of the amplitude of the original wave used 
for the comparisons of theoretical responses. 

Fig. 2 compares the relative performances of the B0BAN3 and fixed methods 
with the traditional method. The RMS errors are taken from within the interior 
of a field of data with highly irregular distribution. Here the standard of 
comparison is the traditional method. The results for the fixed method 
(squares) are superior in the 3S to 8S wavelength range and are inferior in 
the wavelengths greater than 8S relative to the traditional method. There is 
no relative increase in the accuracy of the B0BAN3 analysis (circles) for the 
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short wavelengths because the first pass response for these waves is 
approximately zero. The analyses of the longer wavelengths by B0BAN3 were 
superior to the results from either the fixed or traditional methods. Greatest 
increases over the traditional method were approximately 8 percent for the 4S 
wavelength. 

Comparison of the results from B0BAN3 and the fixed method (Fig. 2) 
reveals the following. Both methods yield similar improvements over the 
traditional method for wavelengths in the range 2S-3S. For AS or greater 
wavelengths, B0BAN3 outperformed the fixed method by 2-3 percent. These 
percentage improvements are for the whole wave and, locally, the improvements 
can be much larger. Now the only difference between B0BAN3 and the fixed 
method is the very smooth first pass of B0BAN3 and the differences between the 
theoretical final responses of the two methods are given by the solid line. It 
is apparent that relative improvements by B0BAN3 in the 2S-5S wavelength range 
can be explained by the first pass of. B0BAN3. Simply add the differences in 
the theoretical response curve (solid line) to the values for the fixed 
method. However, for wavelengths greater than 6S, the relative improvements by 
B0BAN3 can not be explained by the increase in the final response gained by 
including the third pass because the theoretical percentage increases are too 
small. There must exist, therefore, some factor related to the analysis of 
unevenly spaced observations that accounts for the relative improvement of 
B0BAN3 over the fixed method. 

To determine if the factor also accounts for some of the improvement of 
B0BAN3 over the traditional method, we repeated the comparative objective 
analyses for the one dimensional monochromatic waves on a 21 by 21 grid and 
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set the analysis parameters so that analyses would be equally accurate at the 
4S wavelength for B0BAN3 and the traditional method applied to a grid of 
uniformly spaced data. It was necessary to require the equivalence with 
analyses of uniformly spaced data when it was found that the traditional 
method does not filter the short waves to the extent predicted by theory. We 
restored the 4S wave to 75 percent of its original amplitude. The solid line 
in Fig. 3 gives differences between the analyses in percent of the amplitudes 
of the original waves. The traditional method is the standard of comparison. 
Fig. 3 shows that B0BAN3 has the desirable property of restoring more of the 
longer wavelengths (negative percentages) and filtering more of the short 
wavelengths (positive percentages) relative to the traditional method. 
However, B0BAN3 has restored about 2.5 percent more of the 4S wavelength. 

Recall that we set the analysis parameters so that the amplitude response 
of analyses with uniformly spaced data were equivalent. Therefore, the 
differences between the two methods cannot be attributed entirely to 
differences in maximum amplitudes. We therefore define and test another 
measure of the accuracies of an objective analysis - the smoothness of the 
fitted surfaces as one looks normal to the wave front and parallel to the wave 
crests. This is a crucial test because an analysis that better returns the 
amplitudes of the shorter wavelengths as measured by the RMS values for the 
whole wave, as does B0BAN3, may in fact generate more high frequency noise in 
the fitted surface thus giving a "sawtooth" appearance to the analyzed wave. 
It may be better to filter the wave with the traditional method than to retain 
it with a poor analysis especially if there is to be the calculation of 
derivatives of the gridded fields. 

The fitted surfaces should be uniform in the direction of the crests of 
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the one dimensional waves. However, the analysis of irregularly distributed 
data introduces nonuniformities into the fitted surfaces. The slopes parallel 
to the wave crests of the fitted surfaces were calculated for each grid point 
pair located within the grid interior. Since the slopes of the true analytic 
surfaces are zero in this direction, any nonzero values are analysis errors. 
The RMS errors for B0BAN3 were compared with the RMS errors for the 
traditional method (the traditional method was used as the standard of 
comparison) and the results expressed as percent of the amplitude of the 
original wave per grid space. As shown by the dashed line in Fig. 3, these 
percentages are negative for the whole range of wavelengths studied, meaning 
that the slopes parallel to the wave crests as analyzed by B0BAN3 are on the 
whole not as large as the slopes analyzed by the traditional method. 

To provide more detail into the behavior of these two analysis methods, 
we produce in Fig. 4 three pairs of cross sections showing first derivatives 
of the analyzed waves. The first cross section is taken along the crest of a 
wave, the second taken part way between an inflection point and a wave crest, 
and the third taken along an inflection point. The crosses (circles) identify 
the B0BAN3 (traditional) method. There are several ways to interpret these 
cross sections. First, the most accurate method gives the smaller departures 
from zero, the true slope. B0BAN3 almost always gave the smaller departures. 
Second, changes in the sign of the slope identify erroneous extrema in the 
analyzed surface. B0BAN3 tends to produce fewer short wavelength peaks (see 
point 11 Cross Section a). Third, the curvature is a measure of how well an 
analysis method produces smooth surfaces. The greater the "saw tooth" 
appearance of the lines, such as Cross Section b, the larger is the "noise" 
contained in the second derivative. 
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3. Comparisons Between B0BAN3 and Traditional Method Analyses of Real Data. 

It was shown in the last section that B0BAN3 analyses of irregularly 
distributed data are superior to comparative analyses by the traditional 
method in the following ways. (1) B0BAN3 restores more of the amplitudes of 
the desirable long and short waves, (2) B0BAN3 filters more of the amplitudes 
of the undesirable short waves, and (3) B0BAN3 produces analyses with smaller 
interpolation error. But these results are based upon analyses of one 
dimensional monochromatic waves and there remains some question whether the 
improvements are meteorologically significant and justify the added expense of 
a third pass. We therefore perform comparative analyses of meteorological data 
and also compare the first and second derivatives calculated from the analyzed 
fields. 

An intense cylone located over the western United States on 12 GMT, 11 
April 1979 (SESAME I). It was mostly a synoptic event that involved long 
wavelengths and the theoretical response curves predict only small differences 
between the two methods. However, the event also included several jet streaks 
and associated gradients in the height and temperature fields. The data were 
gridded onto a 40x25 mesh with 100 km grid spacing and analyses done for wind 
components, temperature, and height at 100 mb intervals from the surface to 
100 mb. Derivative fields such as divergence, vorticity, and Laplacian of the 
heights were calculated also. Analysis parameters were selected so that the 
final responses for the 2S wavelength were 0.25. 



(13) 


Differences that could be considered as meteorologically significant were 
found mostly in areas of strong height and temperature gradients and also 
within important small scale features such as jet streak maxima and the depth 
of the major cyclone. Comparative results for the 300 mb heights are presented 
in Fig. 5. Gradients were largest at this level and the relative impacts of 
the two methods upon the final analyses were most apparent. Heights analyzed 
by the traditional method and B0BAN3 are shown in Fig. 5a and 5b. The most 
notable difference is that the low over New Mexico (shaded area) has been 
deepened 20 m by B0BAN3, a better fit to the data. Elsewhere, the trough over 
Idaho has been deepened by 5-10 m and the ridge over Wyoming has been 
amplified by 5-10 m. These increases in resolution by B0BAN3 are about 2-4 
percent of the approximate 500 m amplitude of the 300 mb height pattern. 

The Laplacian of the height field is proportional to the geostrophic 
vorticity and is therefore an important meteorological variable. The 
Laplacians of the height fields analyzed by the traditional method and B0BAN3 
(Fig. 5c and 5d) reveal the advantages of B0BAN3 in three ways. First, both 
Laplacian fields show evidence for cyclonic curvature and steep height 
gradients from New Mexico through Nebraska and from Washington through Nevada. 
The amplitude of the pattern over New Mexico is 8-10 percent larger revealing 
that the B0BAN3 analysis has restored more of the desired short wavelengths. 
Second, the patterns calculated from the B0BAN3 Laplacians are smoother and 
show less apparently nonmeteorological short wavelength features. And third, 
the B0BAN3 analyses appear to be less sensitive to boundaries in the data 
field (Achtemeier, 1986). This advantage of B0BAN3 is more tenuous because 
there exist other methods for reducing the impacts of data boundaries (Koch, 
etal. , 1983). 
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In order to identify the short wavelengths in the Laplacian fields, we 
removed most of the longer wavelengths with a high pass filter (Fig. 1) that 
consists of a 2-pass Barnes analysis (fixed method). This filter, applied to 
the gridded data, passes 75 percent of the 2S wave and up to 100 percent of 
the shorter wavelengths. Figs. 5e and 5f show the short wavelengths remaining 
after the application of the filter. The patterns consist of some small 
fraction of meteorologically significant wavelengths and mostly undesirable 
high frequency "noise" waves introduced through the gridding of irregularly 
spaced data. They are mostly 2S waves. An amplitude of 2 units corresponds to 
an error in geostrophic vorticity of approximately 2 * 10 “^ sec - ^. Largest 
amplitudes are found within and around the deep trough. By contrast, Fig. 5f 
shows these waves much suppressed by BOBAN3. Furthermore, the almost 
imperceptible extension of the ridge in the height field over Illinois and a 
weak trough over Nebraska (Fig.. 5a) into a large data void located over 
Missouri has introduced a couplet, a positive and negative center, in the 
Laplacian field (Fig. 5c) over the same area. Consisting mostly of analysis 
error, it has been reduced by approximately 50 percent in the B0BAN3 analysis. 


Excluding from the calculations three rows and columns of grid points 
along the boundary, the variance of the high frequency wavelengths in Figs. 5e 
and 5f was reduced 72 percent by B0BAN3. 

Streamlines calculated by an accurate method proposed by Achtemeier 
(1979) and isotachs of the 300 mb wind (Fig. 6a and 6b) reveal only subtle 
differences. B0BAN3 has increased the jet streak maximum over Kansas by about 
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2 m sec'-*-. Overall, BOBAN3 reduced the RMS differences between the 
observations and vector wind speed interpolated to the locations of the 
observations from the gridded fields by approximately 6 percent. The field of 
differences between the two analyses (Fig. 6c) reveals that B0BAN3 added about 
2 m sec'*- along the full length of the jet streak to bring the wind speeds 
into better agreement with the observations. Elsewhere, the rather large 
differences in wind speed over Idaho and parts of Montana are the result of a 
shift by B0BAN3 of the deformation field southwestward toward the jet streak 
located over Nevada. Since the deformation field was located within an area 
devoid of observations, it was not possible to determine whether the 
southwestward shift by B0BAN3 was an improvement to the analysis. 

Comparisons between fields of derivatives calculated from the components 
of the wind gridded by the two methods (Fig. 7) reveal that B0BAN3 increased 
the magnitudes of the cyclonic shear near the jet streak and also the 
anticyclonic vorticity over Illinois by approximately one unit (Fig. 7a and 
7b). These increases were accomplished without the addition of analysis noise. 
The improvements by B0BAN3 in the divergence (Fig 7c and 7d) appear mostly in 
the divergence patterns. Patterns of divergence calculated from components of 
the vector wind gridded by the traditional method tend to be broken into small 
centers located in data voids midway between observation sites. For example, 
the zone of divergence that extends from Texas into the Great Lakes appears as 
four distinct centers, one each over north Texas, Arkansas, southeastern 
Illinois, and central Wisconsin. Several of these centers are present in the 
B0BAN3 analysis but much reduced magnitudes give the appearance of a smoother, 
more continuous field of divergence. The tendency of B0BAN3 to retain the 
divergence magnitudes and produce more uniform patterns is also apparent in 
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the analysis of the convergence over Iowa and Kansas, and the area of 
divergence that extends from Canada southward through the High Plains to near 
the Colorado border. 

4. Discussion 

Achtemeier's (1987) analysis of the Barnes (1964, 1973) method for 
gridding irregularly spaced weather data found that improvements, particularly 
in the retrieval of short but resolvable wavelengths, were possible if the 
traditional 2-pass method were replaced with a 3-pass method. In this paper, 
the theory for the 3-pass method (B0BAN3) was advanced and tested with 
analyses of analytical and real data. Comparisons between the traditional 
method and B0BAN3 verified the theoretical prediction that B0BAN3 would 
restore some of the desired longer wavelengths, restore more of the desired 
shorter wavelengths, and filter more of the undesired short wavelengths. The 
more specific results of this study are: 

1. Both theoretical results and comparisons between analyses of one 
dimensional analytical waves showed that B0BAN3 would restore 7-9 percent more 
of the desired short wavelengths (3-5S). B0BAN3, upon analysis of the 300 mb 
1200 GMT 11 April 1979 height and wind fields, restored approximately six 
percent more of the magnitudes of short wavelength features such as jet streak 
maxima (2-4 m sec"M and the depth of the major cyclone (20 m). This result is 
in accord with theoretical predictions. 

2. Comparisons between analyses of one dimensional analytical waves by B0BAN3 
and the traditional method revealed that the magnitudes of errors in the first 
derivatives of fields gridded by B0BAN3 were significantly reduced. In the 
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study of the behavior of derivatives of the gridded wind field, it was found 
that the divergences calculated from the B0BAN3 wind field were not 
dismembered into small centers located midway between observations to the 
extent as were the divergences calculated from wind fields gridded by the 
traditional method. 

3. The analysis of the analytical waves also indicated that significant 
improvements by B0BAN3 in the second derivatives were possible. We calculated 
the Laplacians of the heights, an important meteorological variable 
proportional to the geostrophic vorticity. The magnitude of the cyclonic 
center calculated from the height field gridded by B0BAN3 was approximately 
six percent greater than that calculated from heights gridded by the 
traditional method. In addition, B0BAN3 produced a simultaneous reduction of 
the variance of the undesired 2S wavelengths by 72 percent. Given the 
locations of these waves with respect to the locations of the observations, it 
is apparent that they originated mostly through the interpolation of 
irregularly spaced data. The ability of B0BAN3 to reduce short wavelength 
noise in the gridded fields is in accord with the results of the analytical 
studies. 

4. The heights gridded by B0BAN3 were contaminated with less noise of the type 
associated with interpolation across data boundaries. Some of the problems 
inherent in the successive corrections objective analysis in the presence of 
boundaries in the field of data have been described by Achtemeier (1986). 

Whether the modifications of the traditional Barnes successive 
corrections objective analysis into the B0BAN3 method produces 
meteorologically significant improvements to justify the added expense of a 
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third pass depends upon the purpose for the analysis. The traditional Barnes 
method is a highly accurate method for gridding meteorological data. If the 
analysis is to provide a "viewing" of the weather patterns, then the 
traditional method should be satisfactory. If the analysis requirements 
include the calculation of derivative fields and products therefrom, then the 
more precise B0BAN3 method may be preferred. 


As has been noted, meteorological data gridded by the Barnes methods have 
often been used as standards for comparisons with new objective analysis 
techniques. When the objective fields are to be used in higher order 
calculations, such as terms of the differential equations that describe 
atmospheric flow, new objective analysis techniques should be compared against 
B0BAN3 . 
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Figure Captions 


Fig. 1. Theoretical final response curves for B0BAN3 and the traditional 
Barnes method. The responses are identical (0.25) at the reference wavelength 
(2S). Third line: response curve for high-pass filter. 

Fig. 2. (Solid line) Differences between the final responses for the fixed 
method and B0BAN3 for reference wavelength equal to 2S. Relative analysis 
accuracies for the traditional method and the fixed method (squares) and 
traditional method and B0BAN3 (circles) for analyses of irregularly spaced 
data. The relative accuracy is the difference between the RMS errors for each 
analysis pair normalized by the amplitude of the analytical wave and converted 
to percent. 

Fig. 3. (Solid line) Relative analysis accuracies for the traditional method 
and B0BAN3 for reference wavelength equal to 4S. (Dashed line) Relative 
analysis accuracies for derivatives normal to the wave fronts for the 
traditional method and B0BAN3 expressed as percent of the amplitude of the 
original wave per grid space. 

Fig. 4. Cross sections normal to the wave front showing derivatives of 4S 
wave along a) wave crest, b) mid-way between wave crest and inflection point 
and c) inflection point. Crosses (traditional method) and circles (BOBAN3). 


Fig. 5. Comparative results for analyses of the 300 mb heights for 12 GMT, 11 
April 1979 for the traditional method (left panels) and B0BAN3 (right panels). 
Height contours (a and b), Laplacian of heights (c and d), and Laplacian 
residuals after filtering (e and f). Units converted to geostrophic vorticity 
are 10 - ^ sec - * for panels c, d, e, and f. 


Fig. 6. Streamlines and isotachs by a) the traditional method and b) B0BAN3 
of the 300 mb wind, c) Differences in wind speed (m sec - *) between B0BAN3 and 
the traditional method. 


Fig. 7. Comparative of vorticity and divergence for the traditional method 
(left panels) and B0BAN3 (right panels). Vorticity (a and b) and divergence (c 
and d) units are 10 - ^ sec - *. 
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